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ABSTRACT 

We have reduced the data taken with the Spectral and Photometric Imaging Re- 
ceiver (SPIRE) photometer on board the Herschel Space Observatory in the Science 
Demonstration Phase (SDP) of the Herschel Astrophysical Terahertz Large Area Sur- 
vey (H- ATLAS). We describe the data reduction, which poses specific challenges, both 
because of the sheer size of the data, and because only two scans are made for each re- 
gion. We implement effective solutions to process the bolometric timelines into maps, 
and show that correlations among detectors are negligible, and that the photometer 
is stable on time scales up to 250 s. This is longer than the time the telescope takes 
to cross the observed sky region, and it allows us to use naive binning methods for 
an optimal reconstruction of the sky emission. The maps have equal contribution of 
confusion and white instrumental noise, and the latter is estimated to 5.3, 6.4, and 
6.7 mJybeam -1 (1-cr), at 250, 350, and 500 \im, respectively. This pipeline is used to 
reduce other H-ATLAS observations, as they became available, and we discuss how 
it can be used with the optimal map maker implemented in the Herschel Interactive 
Processing Environment (HIPE), to imp rove computational efficiency and stability. 
The SDP dataset is available from |http: / /www.h-atlas.org/| 

Key words: methods: data analysis — techniques: image processing — submillimctcr 
— surveys 
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1 INTRODUCTION 

The Herschel Astro-physical Ter ahertz Large Area Sur- 
vey (H- ATLAS, lEales et al] l2010h is currently the largest 
area open-time key project tha t is being carried o ut by 
the Herschel Space Observatory (jPilbratt et alj f2010). The 
600 hours of observations will map ~ 550 deg 2 of the sky in 
the 5 far-IR and submillimetric bands of t he Photodetec- 
tor A rray Camera fc Spectrometer (P ACS) l|Poglitsch et al.l 
l2010h and SPIRE i|Griffin et alj|2010h photometers. The sur- 
vey of unprecedented sky and spectral coverage, and angular 
resolution, will target sky regions selected in the direction 
of the North and South Galactic Poles, and o f the Galaxy 
And Mass Assembly Survey, the GAM A fields (|Driver et all 
12009(1 . 

The first observations of the project have been com- 
pleted on a 16 deg 2 GAMA target, as part of the initial 
Science Demonstration Phase (SDP). Here we present the 
maps obtained with the SPIRE photometer, and the reduc- 
tion pipeline, which differs in s ignificant asp e cts from th e 
standard pipeline described by iGriffin et alj (|2008l , l201Ch . 
The SDP data are publicly availabl^E and comprise a set 
of three SPIRE photometer maps at 250, 350, and 500 [im 
(discussed in this work), and two P ACS photometer maps 
at 100 and 160 urn (|lbar et alj|2010h . Each map has its cor- 
responding error and coverage maps as separate extensions. 
So urces (point-l i ke and extended) are detected in these m aps 
by iRigbv et~ai1 (|2010l ). and used by ISmith etHl (|2010l ) to 
identify counterparts in catalogues of UV and near-IR se- 
lected sources. 

SPIRE uses arrays of sparsely distributed bolometric 
detectors to sample the sky radiation through scanning at a 
specific angle with respect to the bolometer array axes, and 
the challenge of map making or sky reconstruction is a linear 
algebra problem which has been extensively inves tigated in 
the past (see Section 8 of ICantalupo et al.ll20ld . for a list 
of references) . The bolometric data are combined to provide 
optimal estimates of the sky emission, for which detailed 
knowledge of the detector noise properties is required. Most 
experiments (including the High Frequency Instrument on 
Planck) have only a limited number of detectors, sparsely 
populating their focal planes, but the large format bolo- 
metric cameras, which are now available, pose new chal- 
lenges, as the amount of data is increasing, with potential 
for correlation among detectors, complicating the analysis 
l|Patanchon et alj|2008h . This is of relevance for SPIRE, and 
the detector noise properties need to be investigated in order 
to optimally reduce the data into maps. 

The choice of a map making technique is dictated by 
the noise properties of the data. In this paper we analyse 
the noise in the bolometric timelines, and we use this in- 
formation to generate two sets of maps: one optimised to 
study emission from point sources, and the other to study 
the large-scale structure accessible in the observed sky re- 
gion. The former have been used in the scientific analyses 
of the SDP H- ATLAS data published in the Herschel First 
Results special A&A edition (2010) 

This analysis procedure provides a first investigation 
into the properties of this data-set, and our results can be 
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1 The SDP H-ATLAS data can be accessed from 
http: / / ww w.h-atlas.org| 



The first release of data from the Herschel ATLAS: the SPIRE images 3 



used in future studies to select appropriate map making 
techniques, which may improve on the one used here. How- 
ever, the pipeline we discuss here will provide the necessary 
input to these techniques, in order to optimise their perfor- 
mance. 



2 THE SPIRE PHOTOMETER 

The SPIRE photometer images the sky in three submilli- 
metric bands centred at 250, 350, and 500 um. Each band 
is 30 per cent wide (in A A/ A) and the focal plane is popu- 
lated with a rrays of silicon nit r ide m icromesh "spider web" 
bolometers l|Bock et al. | 1 19961 . 1 19981 ) operated at a cryo- 
genic temperature of < 300 mK. The detectors are cou- 
pled to the spacecraft 3.5-m mirror telescope by smooth 
walled feed horns, and 4.5 K re-imaging optics, which also 
provide shielding from stray r adiation. The linear distance 
of adjacent feed horns is 2/ A (|Rownd et al.| [2003) where / 
is the focal ratio and A is the wavele ngth. Compared to 
its balloon-borne p redecessor, BLAST |Pascale et alj [20091 ; 
iDevlin et al.l [20090 . SPIRE has better angular resolution, 
with point spread functions (PSFs) of 18.1, 25.2, and 36.6" 
full width at half m aximum (FWHM) at 250, 350 and 
500 urn, respectively (| Griffin et al.l [201ol ). Each array has 
respectively 139, 88, and 43 light sensitive bolometers, as 
well as 2 dark bolometers and one resistor. Thermistors on 
each array monitor temperature variations in the substrate, 
due to fluctuations in the thermal bath. The bolometers are 
AC-biased, and the voltage RMS measured across each de- 
tector is proportional to the radiation detected from the sky. 
The constant of proporti onality is the detector 's responsiv- 
ity, which is measured bv lSwinvard et all |2010t ) from obser- 
vations of sources with known spectral energy distributions. 



3 THE DATA 

The SDP H-ATLAS data was obtained in parallel mode, 
with the two Herschel instruments, PACS and SPIRE, op- 
erating simultaneously. The H-ATLAS data set consists of 
two observations - one at the nominal SPIRE scan angle and 
the other with scans at nearly orthogonal direction (~ 85°), 
over the same 16 deg 2 square region centered at coordinates 
RA = 9 h 5 m 30 s , Dec = 0° 5' 0". Each observation is carried 
in ~ 8h by scanning the telescope at the constant angular 
velocity of 60 arcsecs -1 over great circles. The scan lines are 
separated by 155 arcsec in order to achieve good coverage 
for large fields for both SPIRE and PACS maps. During the 
scanning the SPIRE bolometers are sampleqj at 10 Hz. 



4 MAP MAKING METHOD 

Map making is the operation which provides one estimate of 
the sky emission at the observed wavelength, given the data, 
and their noise properties. The data consist of bolometer 
timelines: 

2 For more details please consult SPIRE and PACS-SPIRE Par- 
allel Mode Observers' Manuals available at the Herschel Science 
Centre web pages: http://herschel.esac . esa.int| 



di(t) = A ip (t)S p + m(t) (1) 

where t is the sample's time, i is the detector index in the 
250, 350, or 500 /im array, n,i(t) is its noise, S p is the sky 
emission from the map pixel p, and Ai p (t) is the pointing 
matrix, the elements of which give the weight of the contri- 
bution of pixel p to the detector i at time t. Summation over 
repeated indices is assumed. Throughout this discussion it 
is assumed that S p is the convolution of the sky emission 
with the telescope PSF. 

The well known Generalised Le ast Squares (GLS) esti- 
mator of the sky (e.g. lLuptonlll993l ) is 

S= {A t N- 1 A)- 1 AN~ 1 d (2) 

where all indices have been dropped for clarity, and TV is 
the noise covariance matrix of the data. It is computation- 
ally challenging to apply this equation, especially when there 
are correlations among detectors, and a number of map mak- 
ing techniques have been considered in the literature to deal 
with this problem. However, when the noise is white and un- 
correlated among detectors, this equation reduces to a sim- 
ple, unweighted binning (naive binning), because TV becomes 
a scalar matrix and its eigenvalue is the timeline variance. 

When the timeline noise has a l//-type component, 
high-pass filters can be used to whiten, or suppress a drift- 
ing baseline. This results in an attenuation of the map at 
angular scales larger than A c = /c/k>, where f c is the filter's 
cutoff frequency, and uj the scan velocity. As shown later in 
Section [7] the outstanding stability of the SPIRE photome- 
ter makes it possible to set A c at the map scale. Therefore 
essentially all Fourier modes are preserved in the telescope 
scan direction. Cross-scan modes, however, are suppressed, 
and can only be recovered by scanning at an orthogonal an- 
gle. This can be seen when studying Equati on [2] in Fourier 
space, as discussed in lPatanchon et all (|2008l ). The two SDP 
H-ATLAS observations can be reduced into a map using a 
naive binning, and simple co-addition. Since each map pre- 
serves only one set of modes, the final, coadded map sup- 
presses (by up to a factor 2) the large scales, and therefore it 
is no longer equivalent to the GLS solution of Equation [2] It 
is, however, possible to recover the suppressed modes, since 
all the necessary information is present and the filtering ef- 
fect can be deconvolved. This is demonstrated in Section flOl 
by using the map maker transfer function. 



5 DATA PREPROCESSING 

Each detector timeline is preprocessed using a customised 
version of the SPIRE standard pipeline (|Griffin et al.ll2008l . 
2010), implemented in the Herschel Interactive Processing 
Environment (HIPE), which is run in an interactive mode, 
starting from the Level-0.5 data products (DP). These are 
a set of FITS binary tables grouping together the uncali- 
brated light and dark detector timelines, and all relevant 
housekeeping information. Each table contains data from 
just one scan, so there is a large number of DP tables for 
each observation. 

Pointing information is associated with each bolomet- 
ric sample using the standard pipeline, and the timelines 
are conver ted to flux densi t y per beam by applying the cal- 
ibration of Swinvar d et al.l (|2010h (the calibration is further 
discussed in Section [9}. 
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Figure 1. The tight correlation between the light SPIRE detectors and one of the temperature sensors of each array is shown in scatter 
plots of bolometer timelines vs. temperature in the three top panels. The bottom panels show the same scatter plot, after the temperature 
correlation is removed from the timelines as discussed in Section [6] Detectors in the 250 and 350 /im arrays are stable. The detectors in 
the 500 /im array are less so, and show a residual slow drift which is absent in the dark channel (see Section [6] for a discussion). The 
small gaps in the data correspond to temporal gaps and discontinuities in the timelines. 



Glitches (mostly cosmic ray events, or Bremsstrahlung 
within the instrument) in the bolometer and thermistor data 
are flagged, and replaced with a constrained realization of 
white noise. Two alternative deglitchers are used (alterna- 
tive in the sense that they are available, but not used by 
default in the standard pipeline) to identify these transients. 
One localizes synchronous transients among detectors in a 
given array, and the other is a sigma-kappa deglitcher which 
looks for isolated events. This allows the detection of a larger 
number of low signal-to-noise glitches which could create 
false point-sources in the final maps, if not corrected. 

DC-offset jumps of unknown origin are occasionally 
present in the thermistor and bolometer timelines. These 
jumps cause temporary de-correlation of the thermistor and 
bolometer timelines, so that usual attempts to correct for the 
temperature drift in the bolometer thermal bath (see the 
next section) fail. Unless corrected for, these jumps cause 
large scale-artifacts (stripes) in the final maps, at several 
sigma level. 

At the time of these observations, the standard pipeline 
did not have the capability to search for and correct the DC- 
offset jumps. Therefore transient events have been identified 
by eye and compensated for by shifting the affected timelines 
by an appropriate DC-offset, which is estimated from the 
data. 

The electronic and bolometric instrumental transfer 
functions are deconvolved from the timelines (using the stan- 
dard pipeline) to recover a flat spectral response, and to 
remove filtering phase delays. 

Data products of each scan are linked together to pro- 
vide continuous timelines for each observation (including 
spacecraft turnaround). There is a clear advantage in doing 
this in the analysis which follows, because signal processing 
Fourier methods (filtering, but also power spectral analy- 
ses) perform better, and the unavoidable artifacts in the 



data are confined at the map edges, corresponding to the 
beginning and the end of each observation. Temporal gaps 
in the data acquisition, corresponding to the SPIRE PCAL 
source flashes and other events, are filled with a constrained 
realization of white noise, and flagged accordingly. Flagged 
data are not projected into maps, but filling the gaps with 
noise improves the stability of filters, and reduces pollution 
of good data at the edges of gaps. 

The final preprocessed timeline is a clean (without 
glitches or DC jumps) and calibrated data-set, with associ- 
ated pointing information, which can be binned into a map. 



6 CORRELATED NOISE 

Each light detector timeline contains signal from the sky, 
together with noise. In a large format array the noise has 
two components. One is uncorrelated among detectors, the 
other arises from common signals: 

di(t) =Ai P (t)S p + m(t) + aiC(t) (3) 

where C(t) is a signal common among detectors, 014 is 
a detector-specific constant, and n;(t) is the uncorrelated 
(among detectors) noise, which can have a (time-correlated) 
l//-type spectrum. Here the assumption of only one com- 
mon signal is made, even though in general there could be 
several components with different amplitudes. It is not pos- 
sible to reconstruct sky angular modes with low signal-to- 
noise corresponding to frequencies dominated by l//-type 
noise, because coloured noise does not integrate down with 
time. Since common signals have typically a 1/ /-type power 
spectrum, they need to be estimated and subtracted at 
all sc ales of interest [|Patanchon et all 2008; P ascale et al.l 
12009ft , 

The dominant common signal in the SPIRE arrays cor- 
relates well with the substrate temperature measured by 
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the array thermistors, as shown in the top panels of Fig. [T] 
where representative detector timelines in the 250, 350, and 
500 /mi arrays are plotted against the corresponding ther- 
mistor timelines. The Pearson correlation coefficients be- 
tween bolometers and thermistors are all greater than 99 
per cent. The coefficients a.i are estimated from the slopes 
of these plots. For each bolometer timeline a temperature 
component is obtained from the multiplication of Qi with 
the thermistor. This temperature component is then sim- 
ply subtracted from the bolometer timeline to yield the 
temperature-drift corrected bolometer timeline. To avoid in- 
troducing excess noise, each thermistor is low-pass filtered 
at 20mHz first. Since gradients are not well represented in 
Fourier space, a 5 th order polynomial is fit to the thermistor 
signals and subtracted before filtering, and then added back 
afterwards. 

The corrected timelines are shown in the bottom panels 
of Fig. [T] The 250 /im and 350 fim flat timelines indicate that 
most of the correlated signal has been removed. However, 
this is less true at 500 /im. The temperature-corrected dark 
500 ;Um bolometer (not shown in the plot) has a residual 
similar to the one at 250 and 350 /mi, suggesting that the 
excess signal in the light detectors is optical pick-up. It is 
possible that the 500 fim array is sensitive to the blackbody 
light radiated by the < 4.5 K optics box cavity, which is not 
traced by the thermistor. At 500 um this is 1.5 per cent of 
the peak emission (in vl v of a 4.5 K blackbody), while it is 
negligible at the shorter wavelengths. Although it might be 
possible to decorrelate a second component from the 500 fim 
timelines, this turns out to not be necessary, as shown in the 
following Section. We also notice that this residual is present 
only where the temperature monitored by the thermistor has 
a gradient in excess of ~ 1 nV s~ , measured in the timeline 
units. 

It should be noted that dark bolometers might also be 
used to decorrelate the temperature component from the 
timelines with similar results. This is useful when the ther- 
mometer timelines are not available. However we have not 
addressed in our study the possible issue of light cross-talk 
between light and dark bolometers. 



7 NOISE ANALYSIS 

To understand the noise properties of the clean, 
temperature-corrected bolometers, it is necessary to esti- 
mate noise timelines. This is done by suppressing the as- 
tronomical signal in the timelines, because cirrus and point 
sources have a contribution which can be above the noise 
in individual samples. Naive maps with 6, 10, and 14 arcsec 
pixels, at 250, 350, and 500 /im, respectively, are generated 
by binning the timelines of the two SDP observations. These 
pixel sizes are the default in the standard pipeline. New 
timelines are obtained for each detector by re-scanning the 
map with the telescope pointing solution, effectively simu- 
lating new observations. Since each map pixel contains ~ 10 
samples, the instrumental noise in the re-scanned timelines 
is about 1/3 that of the original timelines. The re-scanned 
timelines are subtracted from the cleaned bolometers to pro- 
vide estimates of detector noise. The noise power spectral 
density (PSD) is shown in Fig. [21 for a representative 250 /im 
detector. 
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Figure 2. The noise power spectral density of a representative 
250 /im detector is shown in the top panel. The solid line is a fit 
of Brown and white noise to the data, showing no evidence of 
steeper power-law components. The noise 'knee' is indicated by 
the vertical dashed line and is in general below 4 mHz for the de- 
tectors in the three SPIRE arrays, which indicates a photometer 
temporal stability > 250 s. The feature at 1.7 mHz (and its har- 
monics) is at the telescope scan frequency. The logarithmically 
binned noise cross-power spectrum between two 250 fim detec- 
tors (diamond symbols) shows that correlations among detectors 
becomes negligible above the noise 'knee'. 



A simple 1// noise model is fit to the PSD 

2 



P(f) 
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where Wo is the white noise level, and /o is the 'knee' of the 
1/f noise. 

The model fit values for fo are below 4 mHz across the 
arrays, corresponding to angular modes larger than 4.2° (i.e. 
the same as the survey angular size), at the scan rate of 60 
arcsec s _1 . The estimates of the white noise level are, u>o ~ 
8.6, 9.3, and 13.2 mjys ' 5 (10 percent estimated variability 
across the arrays) at 250, 350 and 500 um, w hich are com- 
patible with the results of lNguven et alj |2010t) . within their 
quoted uncertainties. 

Residual correlations among detectors are investigated 
by estimating the cross-power spectum between channels in 
an array. An example for two representative 250 /im detec- 
tors is also shown in Fig. [5] It is significant at frequencies 
below the noise 'knee', meaning that the residual 1/f noise 
is highly correlated among detectors. This is not a problem 
for map making because these temporal scales (indicated by 
a vertical dotted line in the Figure) can be safely filtered, as 
they all correspond to unconstrained angular modes (next 
Section). However, we notice that the cross-correlation is 
significantly different from zero at high frequency, indicating 
the presence of broad-band correlated noise of instrumental 
origin. The same level of correlation is observed in the raw 
timelines, indicating that this is not an artifact introduced 
by our analysis. The level of this correlation is estimated at 
1Hz as the ratio of the cross-power spectrum to the PSD, 
which is 2, 10, and 25 per cent, at 250, 350 and 500 /im, 
respectively, and negligible among detectors in different ar- 
rays. The origin of the effect is unclear, however it is not a 
concern for map making because the correlated noise com- 
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Figure 3. The map maker transfer functions in Fourier space are 
plotted for the two SDP scan maps in the left and right panels. 
Fourier modes are attenuated and suppressed in the cross scan 
directions only, where the transfer function is close to zero. Fourier 
modes along the telescope scan directions are not attenuated. 
Combining the two scans allows for optimal reconstruction of all 
the accessible Fourier modes. 

ponent is sub-dominant relative to the uncorrelated compo- 
nent, and because different detectors in an array point at 
different directions in the sky, at any given time. 



8 MAPS 

A 5 th order polynomial is fit and subtracted from each tem- 
perature corrected timeline, to remove residual gradients. A 
line is fit to the first and last 10 sec of data of each timeline, 
and subtracted. Temporal frequencies longer than the 1/f 
'knee' are suppressed by high-pass filtering each timeline at 
the frequency of 4mHz, which corresponds to 4.2° on the 
sky, at the survey scan rate. This preserves all the accessi- 
ble angular modes in one observation (the modes along the 
scan direction) while suppressing the correlated noise which 
would otherwise complicate the analysis or introduce arti- 
facts. The weak filtering adopted does not distort the PSF 
in a significant way, as discussed in the next Section. This 
set of timelines is the final, improved Level- 1 data product, 
which forms the input to the map maker. 

The naive map maker implemented in the SPIRE stan- 
dard pipline is used to bin the improved Level-1 data into 
maps with pixel sizes of 5, 10, and 10 arcsec, at 250, 350, and 
500 nm, with identical astrometry. This pixelization makes 
it trivial to rebin the 250 )im map onto the 350 and 500 fim 
grid, when needed, without loss of information, and Nyquist 
sample the PSF at each wavelength. These pixel sizes are 
smaller than that used by the standard pipeline (6, 10 and 
14 arcsec), and caution should be taken as a few pixels will 
have no bolometer samples (i.e. zero coverage), which means 
they will be left blank. This is not a problem because blank 
pixels are properly flagged in a mask which is also gener- 
ated by the map maker. Three maps are generated: the first 
two are obtained by binning data from the two individual, 
cross- linked observations (Mi and M2), and the other is a 
combined map (M c ) made using timelines from both obser- 
vations. Coverage and noise maps are also generated. 

The two separate SDP scan maps are similar in cover- 



age, and number of samples in each pixel. Therefore, it is 
reasonable to assume similar noise properties in both. By 
subtracting one map from the other, we obtain a jackknife 
which does not contain astronomical signal, and can be used 
to characterise the instrumental component of the noise in 
the maps. A x 2 test of Gaussianity on the jackknife shows 
that the instrumental noise is well described by a Gaus- 
sian distribution, with reduced \ 2 = 1-3, 1-4, and 1.3, and 
Qb? lv) = 15, 7, and 14 per cent, for 250, 350 and 500 /xm, 
respectively (y are the degrees of freedom). 

The instrumental noise RMS, a, is estimated from the 
jackknife pixel histogram, and we find a = 29.6, 31, and 
36 m Jy \J sample (i.e. the standard deviation of each sam- 
ple in a timeline), at 250, 350 and 500 ^m, respectively. 
These values are compatible with the noise estimated from 
the timelines, within ~15 per cent, and are in good agree- 
ment with the noise estimated by the map maker, which is 
provided in the form of a noise map (estimated from the 
variance in a pixel, divided by the number of samples in the 
pixel) . 

We also estimate the confusion noise, the signal arising 
from faint extragalactic sources blended together by the fi- 
nite telescope PSF. Smoothed maps are produced by the 
inverse variance weighted convolution with the telescope 
PSF, and a region of low cirrus emission is selected, from 
which a pixel histogram is formed. This histogram represents 
the combined contributions of instrumental noise, confusion 
noise, and emission from detected sources wit hin the region 
jPatanchon et al.ll2009l ; iMarsden et al.ll2009l ). A Gaussian 
function is fit to the negative part of the histogram only, 
whose width is the square root of the quadrature sum of 
confusion and instrumental noise, the latter estimated from 
the pixel histogram of the smoothed jackknife. 

The 1-cr instrumental noise in the PSF-convolved map is 
estimated to be 4.8, 4.9, and 5.7 mjy beam -1 at 250, 350 and 
500 /im, respectively, with 5 per cent uncertainty. The con- 
fusion noise is, respectively, 5.3, 6.4, and 6.7 mjy beam - , 
with a 7 per cent uncertainty (calibration uncertainties are 
not included) . These can be considered as the relevant values 
when extracting point sources from the maps. 



9 CALIBRATION 

The astrometry has been verified using the point so urces 
detected in the 250 /im SDP maps (|Rigbv et al.ll2010l V We 
produced histograms of the separations in RA and Dec for 
all SDSS DR7 r-band sources within 50 arcsec of the 250 /im 
centroids, in bins of 1 arcsec. We then fit the resulting distri- 
bution of separations using a Gaussian model for the SPIRE 
positional errors and allowing for the effects of ga laxy clus- 
tering in the SDSS data (see ISmith et al.ll2010r ). We find 
a small, < 2 arcsec astrometric shift ( < 1 arcsec uncer- 
tainty), which is compatible with pixelization effects, and is 
accounted for by shifting the maps. A similar analysis using 
stacking of a catalogue of IRAS sources, with positions de- 
termined in the optical and radio, gives comp atible results at 
each w avelength. It is interesting to note that lPilbratt et all 
|2010l) report a 1-cr pointing uncertainty of ~ 2 arcsec, sim- 
ilar to what we find, although stochastic in nature. 

The flux calibration is carried out in the timelines 
as part of the standard SPIRE data reduction (as dis- 
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Figure 4. The 1-D power spectrum of a simulated map of Gaus- 
sian noise with a fc~ 3 power spectrum is shown by the grey solid 
line. The reconstructed map power spectrum (dashed line) shows 
the effect of the timeline filtering in suppressing the larger angular 
scales. Deconvolving this attenuation by using the map transfer 
functions for the two orthogonal scans recovers all the accessible 
angular scales, as shown by the solid black line. The drop off at 
small scales (large k) is caused by the PSF. 

cussed in Section O, and the resulting maps are cali- 
br ated in Jy/beam. We appli ed the flux correctio n factors 
of ISwinvard et all (|201Ch and iGriffin et all (|201Ch required 
by the SPIRE calibration available at the time of this work. 
The 250, 350 and 500 um maps are multiplied by 1.02, 1.05, 
and 0.94, respectively. The map making process effectively 
convolves the sky emission with the pixel window function, 
introducing a gain variation for point sources of up to 5 per 
cent for the 250 /im array. Instead of accounting for this ef- 
fect in the maps, we prefer to include it in the PSFs used to 
retrieve point source fluxes. Gaussian PSFs are fit to maps 
of Neptune, a SPIRE primary calibrator. Th e PSFs are the n 
normalized using the Neptune flux model of lMorenol l|2010l ). 
calculated for the day when Neptune was observed. The 
PSFs, available as part of the SDP release, have FWHM 
of 18.1, 24.8 and 35.2 arcsec, and peak amplitudes of 0.95, 
0.95 and 1.02 at 250, 350 and 500 /jm, respectively. Future 
release of H-ATLAS products will make use of PSFs esti- 
mated from Neptune maps. 



10 SKY RECONSTRUCTION 

In this section we investigate how well the pipeline is able 
to provide an estimate of the sky emission. As the noise 
has already been discussed, it is now necessary to study the 
effect that high-pass filtering the timelines has on the maps. 
This operation is equivalent to convolving the sky with an 
effective filter, or transfer function: 

M„ = T n ®S + noise, (5) 

where n is an index identifying the map (1: first observation, 
2: second observation, and c: combined), T n is the effective 
transfer function, and <S> is the convolution operator. 

3 See http://herschel.esac. esa.int/SDP_wkshops/presentations/IR,/3. 



As mentioned in Section|4l Fourier modes orthogonal to 
the scan direction are lost in the 1// noise, which is filtered 
in the timelines. The two transfer functions for each obser- 
vation are shown in Fig. [3] These have been estimated from 
900 realizations of white noise sky emission, reconstructed 
by the pipeline, and show how the largest Fourier angular 
modes are the most affected. 

The maps preserve a unit response to point sources. To 
verify this, noiseless maps with Gaussian PSFs of nominal 
sizes for the three arrays are generated. The timelines ob- 
tained by re-scanning the maps are then processed by the 
pipeline into new maps. From the comparison between the 
input and output maps, we find that the gain variation is 
below 1 per cent and therefore negligible, compared to other 
sources of error, partic ularly the 15 per cent calibration un- 
certainty estimated bv lSwinvard et all |2010T ). 

The large angular scale Fourier modes which are lost in 
one observation, are constrained by the other. It is therefore 
possible to combine the two maps with no loss of informa- 
tion, when the transfer functions are known. A Maximum 
Likelihood (ML) estimator of the sky is in this case (assum- 
ing a constant and equal noise in each map) 

- _ Mi Tj + M 2 T 2 * 

Ti T* + T% T 2 * ' ( ' 

where the symbols in the equation are the corresponding 
Fourier transforms of the quantities introduced earlier, and 
the complex conjugate is denoted by an asterisk. 

To verify how effective this estimator is, tests were 
performed using a noiseless sky realization of Galactic 
cirrus emission, assuming a k~ 3 power spectrum (e.g. 
iMiville-Deschenes et al.l 120071 . where k is the inverse an- 
gular scale. After reconstructing the emission with our 
pipeline, the input and reconstructed power spectra were 
compared (Fig. 2J) ■ The effect of attenuating the large scales 
is shown by the dashed line. This affects Fourier modes 
which are larger than 20 arcmin on the sky (and therefore 
irrelevant for point sources) . However, the power spectum of 
the ML estimator in Equation [5] (shown by the black line in 
the figure) demonstrates that all accessible Fourier modes 
have been properly reconstructed. 

It should be noted that the scientific analysis of H- 
ATLAS data taken during the SDP phase were constrained 
to point sources, or moderately < 2 arcmin extended sources 
only. Since these sources are not affected by this choice of 
filtering, the naive maps were used, instead of the ML ones 
described here. 

Because detectors can be effectively de-correlated, it 
would be possible to use the well-tested Cosmic Mic rowave 
Background code MADmap (|Cantalupo et al1l2010l ) in the 
future. This is an optimal map maker, which does not cur- 
rently handle detector-detector correlations, and does not 
require an estimate of the map making transfer function 
to compute the GLS sky estimator. The detector noise can 
be made stationary (at least over the time spanned by our 
data-set) and the noise model required by MADmap can 
be efficiently pre-computed and consistently used to reduce 
the whole H-ATLAS ~ 550 deg 2 data-set, when it becomes 
available. A MADmap implementation is now being opti- 
mised in the SPIRE pipeline, and preliminary tests suggest 
_Gtrffihj$KIte.S|)BQffifi&!ptlf reduce the SDP dataset on a single 
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CPU machine, provided that enough memory is available for 
the data. 



11 CONCLUSIONS 

Data from the SDP H-ATLAS observations have been re- 
duced into maps using a customised pipeline. Noise analy- 
sis of individual bolometer timelines demonstrates that the 
SPIRE photometer is very stable, with 1/f residual noise 
dominating frequencies below 4mHz. At the scan velocity 
adopted for this survey, this corresponds to the angular scale 
of the survey region. Naive map making can therefore be 
used to provide the optimal reconstruction of point source 
sky emission. However, to constrain the large scale structure 
it has been necessary to estimate the map transfer function 
by using computationally-intensive Monte Carlo techniques. 

The H-ATLAS maps contain Gaussian instrumental 
noise, in addition to extragalactic confusion noise, and the 
two have roughly the same contribution. Estimates of the 
confusion noise are in agreement with those measured from 
alternative studies. Residual correlations among detectors in 
the same array are found to be sub-dominant compared to 
the uncorrelated (among detectors) noise component, at all 
relevant time scales. 

This map making technique represents an intermedi- 
ate solution until the implementation of MADmap in the 
SPIRE standard pipeline is completely characterised and 
tested. Nevertheless, the improved Level- 1 data products, 
and noise models produced in this work will provide the nec- 
essary inputs in order for MADmap to perform optimally. 
This seams likely to be the method-of-choice to reduce the 
~ 550 deg 2 H-ATLAS survey, as it does not require the ex- 
plicit estimation of the transfer functions, and detectors can 
be efficiently de-correlated at the time and angular scales of 
interest, and moreover the noisy, correlated long temporal 
scales can be filtered without introducing artifacts. 
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